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A. Backgrounä 


In the past several years the Surface Effect Ship (SES) 
has received increased attention in the United States Navy 
and detailed investigations of its sea performance have been 
carried out. The two major categories of SES are the Air 
Cushion Vehicle (ACV) and the Captured Air Bubble Ship 
(CAB). The CAB-type of craft has become of primary interest 
to the United States Navy and two of the constructed 
testcrafts are the Bell Aerospace Systems 100-8 and the XR-3 
of approximately 100 and 3 tons of displacement, 
G sSpecrively. In 1976 the U.S. Navy Surface Effect Ship SES 
100-8 of the CAB-type established a world speed record for 
Mierace type ships of 89.4 knots. Іп April 1976 the SES 
100-B launched a missile at 60 knots. The launch was 


successful and the missile hit its target five miles away. 


In order to investigate the characteristics of the 
CAB-type SES under any design and operating conditions 
without the costly need or modifying the actual craft, a 
Loads and Motions digital computer program was developed by 
Oceanics, Inc. [Ref. 1]. In 1970 the XR-3 testcraft was 
delivered to the Naval Postgraduate School and Leo and 
Boncal converted the 100-B Loads and Motions program to 
represent the XR-3. Since there were substantial design 
differences between these two Ships programming 


nodifications were required in certain subroutines. 





B. Objectives 


The purpose of this thesis is to investigate one aspect 
of safe maneuverapility of these high speed ships, i.e. 
their roll behavior in a turn maneuver. To fulfill this 
objective the sensitivity of the simulated  XR-3 TOLL 
behavior in turn maneuvers as it is represented by the Loads 
and Motions Program is investigated for its dependence on 
changes in sidewall characteristics. The effect of 
modifications in the added mass computation, sidewall 
parameters and force calculations as well as propuision and 
rudder input parameters are investigated. An evaluation of 
the computed performance is obtained by comparison with 


experimentally measured steady state roll conditions. 
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TIF INITIAL RENARKS 


SE PROBLEM OF REPEATABILITY 


It has been this author's experience as well as all 
previous investigators at the Naval Postgraduate School 
using the XR-3 Loads and Motions Program that there have 
been no observed problems related to repeating the results 
СОГ а given input condition. However, during the author's 
first studies of the Loads and Motions Program as listed in 
Menzel's thesis [Ref. 2] an attempt was made to duplicate 
some of the simulation runs given in that reference work. It 
was found necessary as a first step to restore the 
simulaticn program on the IBM 360/67 computer from a card 
source onto a disk.! Using the input conditions stated in 
Ref. 2 it was impossible to obtain identical time 
histories, e.g. for turns. Each time subroutine SIDEWALL was 
read into the computer together with the data input deck, 
thereby  overriding the stored program, an error message in 
the printed output indicated that a division by zero 
occurred in this subroutine. The error message was not 
generated when the stored procran which contained an 
identical version of subroutine SIDEWALL was called. A 
closer look revealed that ln the SIDEWALL-version given in 


Ref. 2 the variable PBAR (plenum pressure) was neither 





1 It should also be noted that results obtained in Ref. 2 
used the IBM 360/67 Release 20.6 system while the results 
obtained in this work used the 21.8B Release with HASP 
installed. 
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defined nor transferred by a COMMON-statement and therefore 
a default value of zero was used during the computation of 
PBHEAD (see Appendix C, SDWL 08880). The missing statement 


РЕДКИЕ РВ = PENNE 


was inserted to subroutine SIDEWALL. With this correction 
the new results were slightly closer to those given in 
Ref. 2. 


The table given below compares output values for roll 
and pitch angle in a turn at 20 kn with constant thrust and 


a 35 degrees rudder step input under calm sea conditions. 


TABLE 1 


— ws s “s a Ña ma sum q s mmf Í ee s me 


without with relative 
PBAR-card PBAR-card Ref. 2 deviation 


Roll angle (deg) 
first peak 7.44 7.39 1 2.6 % 
avrq at 20 sec - 2.65 2.11 d 11.9- 2 


Pitch angle (deg) 
First peak 1.34 153.3 Бү? 9248 % 
avrg at 20 sec 0.99 0.94 928 14.9 7 


Although the differences between columns 2 and 3 in TABLE I 
Might be considered +9 be small as far as magnitudes are 
concerned , the relative deviation is up to 15 % as shown in 
column 4. But what is more important is the fact that the 
roll and pitch angle responses for the uncorrected (plots 1 
and 2 in Appendix A) as well as for the corrected program 
(plots 3 and 4) show an unstable craft behavior for t>40 sec 


Where the pitch angle increases rapidly while the roll angle 


|2 





decreases rapidly approaching the zero degree value. This 
unstable condition is probably due to the 35 degrees rudder 
step input at 20 kn, as it is used in Menzel's study, which 
is physically unreal since it could never be introduced to 
the XR-3 testcraft in an actual run at that speed. 
Therefore rudder angles of up to 15 degrees introduced at a 
rate of 5 deg/sec will be used throughout this work, as it 


is also done on the testcraft. 


Also, regarding plots 1-4 again and not considering 
whether they really reflect the actual craft behavior, it 
should be noted that a short time history ( t<20 sec in 
Ref. 2) could possibly result in wrong conclusions because 
the unstable condition is not evident at that instant of 
time. Therefore, time histories of up to 40 or 50 seconds 


will be shown throughout this study. 
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STEADY STATE CONDITIONS 


Due to the change in subroutine SIDEWALL as mentioned in 


the preceding section the steady state conditions for 


straight runs in calm water have been reestablished Гог 


various speeds and are listed in Table II. 


Taole IL 


amata um "Mh mi eee se qaq 


Steady State Conditions 


DS A AS AO: O O c eee 


first executing the XR-3 
Constant 
keeping the evaluated 


these 


speed 


τος 


Loads 
40 seconds 
enrust 


Steady state values with those given in Ref. 


and  Motions 


Program 


Speed Draft Pitch Plenum Thrust 
angle pressure 
(kn) (in) (deg) (psf) (1 b) 
ο. 6 8.17 1.62 29 οι 400.61 
255 2203 1:11 24.84 349.42 
1320 6.74 0.84 24.84 3935271 
1755 6.41 0.603 24.84 346.45 
200 5.10 0.48 24.84 373.44 
Do 5.37 0.38 24.84 411.53 
гоо 3 59 0.29 24.84 458.62 
2.25 5.48 5 24.84 21331 
HOS 5.28 0.26 24.84 574.43 
These state conditions have been established by 


with 


and then repeating tha run 


constant. À 


2 


comparison 


of 


Shows 


differences especially in the lower and upper speed range of 


again 


anywhere eise since the 


up 


simulation 


program 


Changed since last used by Menzel. 
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to 15 % the cause of which could not be susrected 


been 





tee SIGN CONVENTIONS 


The sign conventions used in the XR-3 Loads and Motions 
Program are identical to those used in Report 71-84 by 
Oceanics Incorporated (Ref. 1j. A rignt handed coordinate 
system is applied to the craft with the positive X-axis 
being measured forward and the lateral Y-axis being measured 
positive to starboard. The vertical Z-axis is measured 
positive downward. Identical signs are used for respective 


forces along those axes. 


The signs of the angles are defined in the following 


manner ο 


pitch angle being positive upward 


(boat noses up) 


roll angle being positive to starboard 


(boat rolls to starboard) 


yaw angle being positive to starboard 


(boat turns to starboard) 


rudder angle being positive to port 
(right rudder, boat turns to port) . 


Zero pitch and roll angle are referenced to the X-Y 


plane parallel to the water surface. 


S 





PENN CAST IGATION OF STIDEWALL EFFECTS 


s muu лл sm s ee mi ee х= À ss х= SS ae SS A AA AA A аа 


ADDED MASS EFFECT 


1. Theory 5 


A basic element in representing the hydrodynamic 
forces by use of slender body theory are the two-dimensional 
sectional values of lateral added mass A22 of the sidewalls, 
which are necessary for the determination of the lateral 
forces, as well as the two-dimensional sectional vertical 
added nass A33 which is also used in determining the roll 
moment. These two-dimensional added mass values are given in 
the ‘Surface Effect Ships  Aero/Hydrodynamics Technology 


Design Manual' (Ref. 3] as 


А22 ο. ΝΑ 


and 


A33 = C *RHO*B*B/8.0 
у 


where E is the lateral added mass coefficient 
C is the vertical added mass coefficient 
у 
RHO is the specific density of sea water 
D is the local draft 


and Beis the Width at the local sidewall waterline. 
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The ES and C values originally selected in Ref. 1 "ere 
у 
those corresponding to high speed with Se апа с=т 
V 
But in order to account for variations in sidewall shape, 


and their influence on the effective lateral added mass, the 
results of the research work described in Ref. 3 led to an 


average value of a: This new coefficient improved the 
agreement between theoretical and experimental data as 


stated in Ref. 3. In the present XR-3 Loads and Motions 


Program a value of e 1s used in accordance with Ref. 1 


L 


and the effect of letting Wasa will be investigated in 


later parts of this thesis. 


Considering the equation for the vertical added mass 
A33, a major difference has been found between that one 
given in Ref. 3 and the one used in the XR-3 Loads and 
Nocsons Program (e.g. Ref. 2), where in accordance with 
Ref. 1 (1971) there was an addition2l PI-factor in the A33s 
equation. The letter 's' indicates that the computation of 
the added mass is done at the craft's stern. Since Ref. 3 
(1976) reflects the experimental and theoretical work being 
done based on Ref. 1, the A33s equations in subroutines 
SIDEWALL and SIDETAB have now been changed to the new form 


given above. 
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The change in the simulated XR-3 performance due to 
the reformulated  A33s 
TABLE III 
Straight runs in calm sea for 
νυν. and 0.8) 


conputation can be observed from 


giving the new steady state conditions for 


various speeds which are 


identical for lateral added mass 


Ececfficients. 


TABLE III 
Steady state conditions 

(new_A33s) 

Speed Draft Pitch Plenum тет 
angle pressure 
(Xn) (in) (459) (ps£) (15) 
150 6.80 0.86 24.84 3230230 
17.25 6.16 0765 24.84 SUD 
20.0 беш] ooo οι. eu SS 
2:299 9.9/2 937 24.84 ΠΠ 552 
25.50 5571 0.30 24.84 соо 
2 5 5.53 νυν οὐ 24.84 555) 
30.0 5.40 0.28 24.84 515206 
Comparing these steady state values with those previously 


given in TABLE II in Sect. II.B one finds in general for all 


speeds that 


агас лаз 1псгеазеа by about 0.05 in (0.8 %) 


- pitcn angle has increased by 0.02 deg (2 to 7 %) 


- thrust nas slightly increased by less than 0.5% 


18 
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J. Effect in a Turn maneuver 


The effect of the change in the formulation of the 
vertical added mass in a simulated turn maneuver has been 
studied next. After a 5 sec straight run at 20 kn in calm 
Sea a port rudder deflection was introduced at a rate of 
5 deg/sec and then kept fixed at 15 degrees resulting in a 
euren to port. The deadrise forces are computed at the 
transom. In TABLE IV are shown the steady state conditions 
and the roll moments contributed from the bow seal (FKBS), 
stern seal (FKSS), sidewalls (FKSW), rudder Cee 
propuision system (FKP), aerodynamics  (FKAED) and plenum 
pressure (ABPB*PHI*Z). 


old A33s new A33s deviation 

pitch angle (deg) 0.48 0232 837m 

roll angle (deg) 230 1.98 - 1.0 % 

Graft (in) 5.86 5.94 1.4 % 

speed (Kn) 19225 19.27 О 
moments: 

FKBS -391.1 -386.3 - 1.2 % 

EKSS 556 - 3.6 0.0 4 

FKSW 191.0 20220 эү > 

FKRud = 97.5 OO -12.8 % 

FKP 02 022 0.0 & 

FKAed - 54.4 = 539 1.7 % 

ABPB*PHI*Z 3595.4 sod. 9 {ο 
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From TABLE IV it can be seen that, due to the 
modified A33s computation, the counteracting roll  mcments 
from the sidewall (FKSW) and rudder (FKRUD) changed by about 
6 and 13 €, respectively, resulting in only a little effect 
in generating a tendency toward a smaller roll angle which 
is directed out of the turn. Other major contributions to 
the roll angle are from the bow seal (FKBS directed into the 
turn) and from the plenum pressure (ABPB*PHI*Z) directed out 
EU the turn. The corresponding plots (plots 5 ~ 8 for old 
ÀA33s, plots 9 - 12 for new A33s) show that with the new 
equation for A33s the pitch and roll angle response curves 
are less damped than they were uSing the former equation for 


vertical added mass calculation as would be expected. 


It should be noted that the roll moment due to 
rudder is fairly large and changes quite Significantly while 
the moment contributed by the propulsion is very smail and 
constant. The cause of this will be investigated in 


ENapter IV. 
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B. DEADRISE ANGLE 


1. Background 


= mio mM M  —— = 2 = 


SIrEgeerErrectsseships. (SES) of the captured air 
bubble type (CAB) have rigid sidewalls immersed into the 
water, thereby - together with the bow and stern seal - 
capturing the air in the plenum chamber and reducing air 
leakage. Since the sidewalls of the XR-3 are not uniform in 
cross-sectional shape throughout their length but are curved 
on the outboard side near the bow similar to a single hull 
ship, the expression "sidehulls" would be more appropriate, 


КЕШ сап be found in modern literature. 


The importance of the correct understanding of the 
effect of the sidewalls on the craft's performance is 
readily seen from a report on SES Research with the XR-1 
testcraft (Ref. 5]. As reported, the first test series were 
carried out successfully with a 45 degrees sidewall deadrise 
angle, which is the angle between the horizontal plane and 
the ship!s outer sidewall surface. In 1964, after a 
modification o£ the craft to 60 degrees deadrise angle and 
articulating seals, an unstable roll condition occurred 
during a turn at 35-37 knots, a maneuver performed many 
times before. The testcraft heeled out of the turn , nosed 
down due to the retraction of the bow seal, continued its 
Seward roll motion aad then flipped over. The report, 
however, does not mention the rudder action generating the 
Burn. After this accident, the deadrise angle had been 
changed back to 45 degrees. From this experience the 


surface effect ship's roll stability should be expected to 
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be sensitive to the deadrise angle and an investigation of 


this sensitivity is reported in the following sections. 


2. Forces due to Deadrise Angle 


The deadrise angle, as pointed out by Ref. 3 and 4, 
has a major effect on the craft's roll motion and therefore 
is an important design consideration. During a turn the 
craft is desired to heel inboard, thereby minimizing the 
apparent transverse acceleration (a coordinated turn). For 
typical SES designs, it is not possible to achieve perfect 
coordination, but it is possible to achieve nearly flat 
mit~ns if the craft is designed properly. The principal 


factors affecting roll stability at speed are 


* Sidewall geometry 


= сеа characteristics 


* vertical location of the center of gravity. 


The forces acting at the center of gravity station of the 
starboard sidđewall in a port turn are shown in Figure 1. 
The centrifugal force generated acts away from the center of 
the turn. This force must be counteracted by hydrodynamic 
forces which are created by yawing the ship into the turn. 
To an observer aboard the craft, it would appear as if che 
craft is sideslipping. Wave buildup on the outboard sides of 
the sidewalls increases the pressure there, while only small 
pressure changes occur on the inboard side of the opposite 
Sidewall generating a small (negligible) vertical force. 
The principal force arises from the outboard side of the 
sidewall. The direction of this force depends on the slope 
of the deadrise angle. If the resultant force passes above 


the center o£ gravity (solid line in Fig. 1) a restoring 
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moment results, tending to heel the craft into the turn 
thereby improving its roll stability. If the deadrise angle 
is chosen to be larger (e.g. 70 deg, shown dashed), the 
resultant force passes below the center ОГ gravity. 
resulting in a moment tending to heel the craft out of the 


Burn. 


Thus, once the SES geometry and the exact vertical 
location of the center of gravity are known, an approximate 


Check for stability in a turn can be made. 


= 


Iure di - GENERAL ROLL STABILITY CHECK IN A PORT TURN 
Ra ERTEIGBE FORCE; C.G. CENTER OF GRAVITY) 
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From the construction data entered into the sidewall 
table (subroutine SIDETAB) Figure 2 has been drawn showing 
the cross sections at several selected stations of the 
Starboard sidewall. Since in most simulation runs (which 
will be discussed later) draft was Less than 7 inches, the 
deadrise angles for all stations immersing into the water 
have been calculated by straight line approximation and are 
Seven in TABLE V. The statiors are counted from bow 
(station 0) to stern (station 28) and are not equally spaced 
for station numbers less than 11. The center of gravity has 
been located experinentally at station 14. The average 


deadrise angle of all stations listed is 64.3 degrees. 


TABLE V 
Deadrise angies (deg) for ali stations 
Straight line approximation for 7'' draft 
Station Deadrise angle Station Deadriss angle 

5 67.6 17 62.1 
6 70.8 18 64.0 
7 57.9 19 65.3 
8 54.4 20 66.7 
9 53.7 21 68.1 
10 53.7 22 6.96 
11 53.7 23 70.3 
12 55.6 2% Tale 
13 57.3 25 nou 
14 58.5 26 73.8 
15 59.7 27 76.9 
16 60.9 28 78.8 
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4. Effect of Deadrise Angle 


From the previous discussion the surface effect 
ship's stability behavior is expected to be sensitive to tne 
sidewalls’!  deadrise angle. In the XR-3 Loads and Motions 
Program listed in Ref. 2 the deadrise angle at the  transom 
(78.8 deg) had been used for the computation of the deadrise 
rorces. This calculation was implemented into subroutine 


SIDEWALL using the following statements: 


FYH (J)z-A22S*U* (V*XSS*R-ZS*P) 

° 

° 

CTNDR=0.0 

IF(DSS.LE.0.0) GO TO 22 

CTNDR=(ES-38(1)) /DSS 

IF(THETA.LT.0.0) CTNDR=0.39391 
22 CONTINUE 

FZHOLD (J) =FZH (J) 

FZHDRP (J) =PM1*FYH (J) *CTNDR*PROMO1 

FZE (J) =FZH (J) +FZHDRP (J) 

° 


FK= (FZH (2)-FZH (1) )*YSW+ FKD-FY*ZS 


The projection of the lateral force FYH(J) at each sidewall 
(J=1 or στα I2 τος Starboard side) produced the 
respective ccmponent of vertical force FZH, where FZHDRP is 
the projected force and CTNDR is the cotangent of the 
deadrise\angle. Tke roll angle had been already taken into 
ιτ in the computation of DSS (draft at the stern). The 
factor PRCMO1 simply provided a means to arbitrarily change 
the vertical projectec force for studies undertaken in 


Ref. 2. The factor PM1 introduced the necessary sign change 
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in the sidewall force computation being dependent on the 
craft's side (EM1=-1 for port, PM1=+1 for starboard side). 
The roll moment FK was partly determined from the vertical 


forces of both sidewalls. 


Comparing the deadrise angles given in TABLE V it is 
seen that the transom deadrise angle of 78.8 degrees is 
unique and the largest along each sidewall and no reason 
could be found why just the transom deadrise angle should be 
used for the vertical force computation. This fact rather 
came up when Leo and Boncal (Ref. 8) scaled down the 
Брига сл оп ргодгап оғ the 100-B testcraft to create the Xa-3 
Loads and Motions Program. poz ene 10055 testoratt che 
Sidewall cross sections are uniformly shaped throughout most 
Meats length. This uniformity does net occur for the XR-3, 
as shown in Figure 2. Therefore it has been supposed that 
another deadrise angle which is more representative for all 
angles existing along the XR-3 sidewalls, e.g. the deadrise 
mabe at the center of gravity Location (station 14), could 
possibly be more appropriate for the vertical force 
memputation. To investigate the effect of this supposition, 
turn maneuvers in calm sea at 20 knots and various rudder 
angles have keen simulated for both deadrise angles and 


lateral added mass coefficients. 


TABLE VI 


a sm Àm Qam Q чы = 
ED aut “mamam το αν щщ ci uma Qamqa; Q= с сз s = = πα. == ο . 


E “mm jm Í s O A w  — = — = =— — = w 


Rudder Deadrise force computation at Experimental 
angle transom center of gravity tes Hen Le 
=0. C z0.4 C =0.8 data (Ref.d) 
h h 
5 DSZ ο... DIIS ου. ο 
9 15:0 0.78 (EL ο 28 
12 1250 Te ES AA ο. 55 
15 Ss Tesa unstable 1.50 
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Table VI shows the steady state values for the roll 
angle which are compared with those measured experimentally 
in 1974 and documented in Ref. 6. Since the XR-3 Loads and 
Motions Program under investigation represents the XR-3 
craft configuration as it existed in 1974, Ref. 6 may serve 
to check whether the results produced by certain changes in 
the simulation program are favorable or not. It is not 
anticipated to exactly match the results obtained from the 
simulated runs to the measured values since for the 
testcraft data the following precisions for the measuring 
Mevrces are stated in Ref. 6 : pitch aná roll angle +0.50, 
rudder angle +1.0° and speed +1 kn. From TABLE VI it is 


obvious that the steady state values for roll angle using 
E^ as the lateral added mass coefficient and the center 


of gravity deadrise angle for the vertical force computation 
are closer by 25 to 40 % to the measured angles than are the 
corresponding roll angles considering the transom deadrise 


angle. 

Using the center of gravity deadrise angle and the 
lateral added mass coefficient m? as suggested in 
Ref. 3, the agreement in steady state roll angles could be 


improved by another 3 to 10% for rudder angles up to 
12 degrees. This improvement has also been reported in 
Ref. 3. Comparing the roll angle responses for turns 


generated by a 12 degrees rudder angle to port with cau 
(plots 21 and 22) and αν (@lots 23 and 24) it “is round 
that the smaller lateral added mass coefficient generates a 


little higher peak roll angle (1.479) than does the larger 
севе (1.279). Both responses show the same number of 
oscillaticns (eleven cycles) until they die out. The larger 


lateral added mass coefficient generates a small negative 
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\ 


roll angle shortly after the port rudder motion has been 
introduced. However, for a 15 degrees port rudder angle the 
Simulaticn program showed an unstable response with 
increasing amplitudes of oscillation for pitch and roll 
angle (see plots 17 and 18). From this the lateral added 
mass coefficient being 0.8 does not seem to be appropriate 
for the XR-3 Simulation if both sidewalls are considered for 
the vertical deadrise force computation. The effect of 
using both lateral add2d nass coefficients will be 


Sonsidered again in Sections III.D and E. 


From TABLE VII it can be seen that the roll moments 
contributed by the bow seal (FKBS), plenum pressure 
(ABPB*PHI*Z), rudder (FKRUD) and sidewalls (FKSW) change the 
most. Comparing the effect due to the change from transom to 


center of gravity deadrise angle it is found that 

- bow seal inward effect decreases 

- plenum pressure outward effect decreases 

- rudder inward effect decreases 

- sidewall outward effect increases 

- outward roll angle decreases. 

The time histories for simulated turn maneuvers with 
15 degrees port rudder angle using center of gravity and 
transom deadrise angle (TABLE VII) are shown in Appendix A 
ο ροτς 13-16 and 9-12, respectively. Comparing these 


Pees the effect of the change from transom to center of 


gravity deadrise angle and us can be observed as 
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- pitch angle response being more damped 


- roll angle response being more damped 


- pitch and roll rate responses being more damped 


and having smaller peak values. 


Steady State Values for selected Deadrise Angles 


A ma aras < = O En SS m s m RS En me SEE Ten Q= mm m s mm s s eee = e 


Deadrise Angle Transom ES OL Ges arbitrarily chosen 

(degrees) 78.8 5979 316 28.7 
Pitch angle (deg) 0452 0.45 059/2 0.6ü** 
Roll angie (deg) 1-38 oe s ο ο ΟΞ 55Ξ 
Draft (in) 5.94 5:51 6 07 bea 
Speed (kn) 19227 19.37 quor 1827 
FYSW =824 0 - 21.9 ugs 207,0 
FYRUD 40.0 ЗР η D O 
FYAED Ошо το. Saco ed 5 
R*V*AM 80325 80922 0319 8 πο, Ὁ 
BKES ΘΠ 2336.4 - 68.8 71.6 
BASS 5 ΞΘ“ «EOS 025 
FKSW 29220 Zu 2 21 17229 184.0 
FKRUD πο. 305009 pu - 53.0 
FKAED 5-5 -. 55... 9 ο ο. = 46.2 
FKP ΙΕ 0.14 0.03 AOS 
meee Pol *Z 3519 2155 58.6 = эу? 
ПОН (Р/5) ОЕТ τ. SOLIDE =107/-88 
BYD SAS -651.4 -637.3 -607.0 
Note * = still decreasing 

** = still increasing, not quite steady state 
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The reason for the sidewall roll moment and thereby 
the net roll effect not changing more pronounced (TABLE VII) 
is that not cnly the lateral force FYH must be considered in 
Order to determine the force F effecting roll stability 
КЕСЕ, ТТТ.В.2, Fig. 1) but also the lateral drag force FYD 
has to be taken into account. Considering the XR-3 geometry 
using the center of gravity deadrise angle, as it is shown 
in solid lines in Fig. 3, the lateral force FYH can be seen 
to generate the vertical projected force FZHDRP and the 
force F, which is directed well above the center of gravity. 
But since FYD has to be added vectorially to the force F and 
its strength being several times larger than FYH (TABLE VII) 
the new resultant force Е! determining the craft's roll 
behavior is directed well below the center of gravity 
(although F being directed above it). Thereby an outward 


roll angle is generated. 





444 289 ὐ 5. 


Figure 3 - DEADRISE AND DRAG FORCES IN A PORT TURN 
(CG= CENTER Or GRAVITY, WL= WATER LINE) 
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Also from the 


center of gravity versus the 


mig. 3 


force generating the inward or 
The of 


corresponding forces are shown in 


seen. center gravity 


effect 


transon 


deadrise 


of 


deadrise angle on the 


considering the 


outward roll moment can be 


angle and its 


solid or dashed lines with 


of The 
lines, generates a force F 
ОЕ 

gets 
the 


F being directed above the center Gravity. transom 


deadrise angle, shown in dotted 


being directed below the center 
Toll 


pectorially FYD. 


απασαινσ producing 
by 


Stability 


an 
which reinforced 


this 


outward noment adding 


From general Check 


yields a smaller outward roll moment and angle using the 


center of gravity versus the transom deadrise angle. This is 
in agreement with the results given in TABLZ VI. 
deadrise 


Next, some Simulation runs with 


chosen 


атр ават у 
angles such that small positive and negative steady 
State roll angles result are presented in order 
find 


necessary to simulate flat turns, 


ПОЕ 55 


out which deadrise angle would generate the moments 
the 
trend of the change in the respective moments which has been 
эг gravity 
The value of the lateral added 


EET coff:cient in these runs was 0.4 as it was in all runs 


and, “second, verify 


observed for the change from transom to center 


deadrise angle in TABLE VII. 


listed in TABLE VII. Considering the roll noments over the 
range of deadrise angles from 78.8 to 28.7 degrees, tne 
moments due to the bow seal, stern seal, propulsion and 


plenum pressure change consistantly in nagnitude while those 
due to sidewalls, rudder and aerodynamics reach sone extreme 
value and then increase again. state 
τοῦ] listed 


arbitrary chosen 


Regarding the 
νους VII; 


deadrise 


steady 
it can be seen that the 
το ра 


that would generate zero 


angles en 


angles happen neariy 


synmetric about the deadrise angle 


nC angle (flat turn) under these simulation conditions. 


Keeping this symmetry in mind one finds that the values of 
the roll moments due to pow seal (FKBS), stern seal (FKSS), 
propulsion (FKP) and plenum pressure (ABPB*PHI*Z) are also 
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Symmetric with respect tS zero moments. For propulsion and 
plenum pressure the sign of the moments are identical to 
that of the roll angle while for bow and stern seal moments 
the signs are reverse. From the results given in TABLE VII 
with positive steady state roll angles Figure 4 has been 
drawn where the dotted line represents the craft's initial 
state, the dashed line represents the craft's steady state 
considering the transom deadrise angle and the solid line 
represents the craft's steady state (smaller roll angle) 
considering the center of gravity deadrise angle for the 
Aoi cal force calculation. Same markings appiy to the 
moment vectors shown from which their change in magnitude 


for both deadrise angles can be seen. 





__-Rup 
ed 


Figure 4 - ROLL MOMENTS IN A PORT TURN 
(ROLL MOMENTS AND ANGLE NOT TO SCALE) 


Legend: AED Aerodynanic CG Center of gravity 
ES Bowseal PP Plenum pressure 
SS Sternseal SW Sidewalls 
RUD Rudder P Propulsion 
Deadrise angle: - - - Transon 


Center of gravity 
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Since the simulation run with the center of gravity 
deadrise angle being used for the vertical force computation 
gave results that were more favorable (better damped time 
histories and steady state values closer to those measured 
experimentally) this change in vertical force computation 
was implemented into the XR-3 Loads and Motion Program. 
Starting with the deadrise angle given in Table Y for draft 
up to seven inches, again straight line approximations for 
the effective deadrise angle at the center of gravity 
location (station 14) for larger drafts have been performed 
NSSulting in an almost linear relationship. This part of 
subroutine SIDEWALL is now valid for any draft and contains 
the following statements which have been used in this form 


E)! simulation runs to follow in this study: 


DDRANG=0.0 

IF (DS.GT.0.5833) DDRANG=(DS-0.5333) «0.0629 
DRANG=1.021+DDRANG-PMI*PHI 

CTNDR=COTAN (DRANG) 

FZHOLD (J) =FZH (J) 

FZHDRP (J) =PM1¥*FYH (J) FCTNDR*PROMO1 

FZH (J) =FZH(J) +? ZHDRP (J) 

e 


FK= (FZH (2) -FZH (1) )*YSW+FKD-FY*ZS 


It should be noted that vertical projected forces 
from both sidewalls are used in computing the roll moments 
in the XR-3 Loads and Motions Program. Section III.E of this 
thesis will be concerned with the effect of changing the 
simulaticn program to consider the deadrise projected force 
only from the outward sidewall in a turn maneuver as it 15 


HENScussed in Sect. III.B.2 . 
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EXC CROSS-FLOW DRAG 


mathie preceding section the lateral force FYD which is 
due to cross-flow drag has been found to be rather large 
compared to FYH, thereby changing the direction of the 
resultant force significantly. FYD is computed in 
ENuDroOoutQine SIDEWALL by suanming the contributions DF of all 


stations: 


BEST, I) = -HRHO*CDSW*VREL*ASS(VRDEL)*DSWAV(L,dJ) 


where DF (I,J) is the lateral drag contributed by 
the j-th aleman of each (i-th) sidewall 
НЕБО is РБӨ/2. 
VREL is relative velocity 
DELX is the incremental distance used for the 
исе ра: calculation along each sidewall 
DSWAV is the average wetted draft over one 
incremental distance 
Slswers tiewcnoss-tlow drag coefficient for one 


sidewall . 


The LOSS flow drag cCoerficient being used in the present 
XR-3 Loads and Motions Program is the value corresponding to 
mee for flew normal to a long flat plate being CDSW=1.28 . 
BM points out that an investigation of cross-flow drag 
forces and coefficients (before the vertical added nass A33 
had been reformulated as discussed in Sect. III.A.1) showed 
better agreenent between tneoretical and experinental 
results if CDSW=2.0, representing the normal агаа 
coefficient for a sharp flat plate, had been used instead. 
After the reformulation of A33 the agreement between 


modified theory and experiment had very much improved, so 
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that the former drag coefficient could be used again. Since 
in this thesis work the change in A33s did not effect the 
roll behavior significantly (TABLE IV), some turn maneuvers 
at 20 kn and 15 degrees rudder angle nave been simulated for 
various drag coefficients in order to investigate which CDSW 
would better match simulated and experimental steady state 
roll angle. Prom TABLE VII a decrease in roll angle and 
roll moment is desireable and therefore the cross-flow drag 
force FYD should be decreased via reducing CDSW. The steady 
state values are shown in TABLE VIII with the lateral added 


mass coefficient being unchanged с =0.4. 
h 


ας шшш самы; cee шо 


CDSH 1283 12.16 1.0 Ў 
Pitch angle 0.43 0.43 0.44 0.46 
Roll angle 15955 1.48 1738 Tag 
Draft J18 7.00 5.55 Ss 91 
Speed ¡LS πο 19.31 wa 
EKBS -343.0 US -308. 1 22909 
EKSS 27 O Zee ο 0 
FKSW 20728 1750 DTS 40.0 
FKRUD = - 47,9 5-5 80.4 
FKAED = 56 26 ЕЯ 285146 - 67.6 
FKP 0.14 0.14 oq 03] 
ABPB*PHI*Z 276.1 264.6 246.5 201.9 
ες (P, S) -4 3/ -125 -45/-125 -48/-124u -57/-121 
PYD -654.9 5 19.5 -561.0 -436.0 


Comparison with Ref. 6 (Table VI) shows close agreenent 
INTO] angle (1.360) for the run with CDSW=1.0 . But these 
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trial runs have been executed with arbitrary smaller drag 
coefficients because an increase in drag coefficient (e.g. 
CDSW=2.0 for sharp flat plate as in Ref. 3) would have 
enlarged the disagreenent between simulated and experinental 
EN Also, from plot 29 for CDSW=1.28 and plot 31 for 
CDSW-1.16 it is seen that a smaller drag coefficient results 
in more damping in the roll response, in this case six 
versus nine cycles of distinguishable oscillations. Since 
the choice of CDSW=1.0 cannot be declared to be appropriate 
for the actual shape of the XR-3 sidewalls it will not be 
considered any further in this study.  CDSW-1.16 generated a 
4 % change (about 15 % is desired) toward the experimentally 
measured roll angle and its effect will be investigated 
together with the present cross-flow drag coefficient 
CDSW=1.28 in the studies to follow. 


Shown below in Figure 5 is the dependence of steady 
state roll angle in turns generated by a 15 degrees rudder 
Saale to port at 20 kn on the cross-flow drag coefficients 
used in the simulation runs (TABLE VIII). 


Roll angle 
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O THRUST MAPPING 


In turn maneuvers as they have been discussed in the 
preceding sections another important point to be considered 
is the change in thrust on both engines. Since the engines 
mounted on the XR-3 provide constant power for a given 
ито е setting, the actual thrust delivered during a turn 
Maneuver may vary even though the throttle remains fixed. 
In all simulation runs shown in this thesis until now 
constant thrust has been used. In Bef. 6 experimentally 
measured reduction in thrust is given for turn maneuvers at 
20 kn and various rudder angles. From this and Ref. 7 
providing more detailed information the thrust-rudder 


dependence was as listed below. 


TABLE IZ 


ss Qw mam q Qe 
нч ew ew ee ce eee ee ee ee - х= ee E чы CO CA сш: es E S =s s s sn - = ae SS a 


ar <p ma A O A O A O ce шшш шшш — = = сл 


Rudder angle corali reduction in thrust (Ret, 6) 
(deg) (2) 
220509 
zu, 
2 >= 4.3 


Γ.Ε τί has been used for thrust mapping (Block 16 ος 
the data input) in the simulation runs. Due to 
unavailability of more precise data it has been assumeä that 
F total reduction in thrust is contributed from both 
engines in equal amounts. 3ut this may not be quite true, 
since due to the craft's roll angle a difference in 


immersion depth for both engines could result in unequal 
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changes in thrust. TABLE X shows the steady state values of 
these simulations considering two cross-flow drag 
coefficients, deadrise force contributions fron both 
Sidewalls computed at the center of gravity as well as both 


lateral added mass coefficients gael? s 


TABLE_X 


mn ee eee wp waar a “a “W - Se ee ce ee eee se cee ee eee = =. = = 


Ine20 kn turns at various rudder angles 


Rudder angle Roli angle (deg) Speed (kn) 
(deg) simulation Ref. 6 Simulation Ref. 6 
CDSW=1.28 
ОС 010725) 0.09 23009. 2) 1953 
ОСО е 0028 T2- 919720) 13255 
12 099006597) Us 15, 515-50) 1922 
CDSW=1. 16 
on» e. 23) Ue, TF2 (192) OR 
086710557) 028 ο ο ας ο) 1955 
12 0.94(0. 84) 0.58 | 11555) 192 


Comparing the steady state roll angles obtained if 
thrust mapping is used with those listed in TABLE VI (no 
thrust mapping, deadrise force computation at the center of 


ΕΥ), a change in roll angle by 6 to 14 % for ul and 

by 14 to 20 % for C =0.8 towarä the measured angle can be 
h 

observed. The steady state roll angles found in the 


Sega tion runs still differ from the roll angles given in 
Memeo by a factor of about two and three but can be 
expected to be a little cioser to these if proper thrust 
mapping, i.e. different maps for port and starboard engine, 


is used. By this also the steady state velocities which 
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differ by about 2.5 % can be expected to come in tetter 


agreement, i.e. to increase for smaller roll angles. 


Regarding the roll angle time responses for sinulation 
runs uSing both lateral added mass coefficients (0.4 and 
0.8) and 12 degrees rud3er angle the Following 
characteristics have been found: the runs using thrust 
mapping (Table X, plot 33(35)) show almost the same 
transient behavior as do the runs without thrust mapping 
K D] VI, plot 21(23)) with about ten cycles of oscillation 
before the transients die out. But the runs using thrust 
mapping need quite a long time (t>45 sec) before they reach 


the steady state roll angle. 


The effect of using a larger lateral added nass 
coefficient (0.8) can be observed in plot 35. Comparing 
miso) with plot 33, the roli angle response in plot 35 shows 
a smaller peak roll angle (1.259 versus 1.470) and the 
transients have smaller amplitudes. But both die out after 
ten cycles. This efiect has alseady been noted in Section 
INL P.14. 


Investigating the effect of changing the coefficient 
CDSW=1.28 to 1.16 by comparison of the roll responses shown 
EN ots 33 and 37 or 35 and 39, a significant damping 
effect as has been found previously in Section III.C can not 


be observed here. 
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Pee ODEADRISS FORCE OF OUTWARD SIDEWALL 


ies Dart Of the investigation on how to effect the roll 
behavior of the testcraft simulation is again concerned with 


the deadrise angle. 


In the present XR-3 Loads and Motions Program deadrise 
forces for both sidewalis are computed with the deadrise 
angle appropriately corrected by the roll angle for port and 
Starboard sidewalls. The idea benind this, e.g. in a port 
Eun) xs to consider a pressure buildup (large РҮН (5)) on 
the outboard side of the starboard sidewall and a reduction 
in pressure on the outboard side of the port sidewall (small 
ШЕН (р) у. Thereby an upward vertical force FZHDRP (with 
negative sign) is obtained for starboard side, while a 
Hownward vertical force (positive sign) is computed for port 
side. Both vertical forces are then added to the buoyancy 
ο Of their side. The forces obtained by this approach 


are shown in Figure 6. 


In Ref. 4, however, a different philosophy in regarding 
the acting forces is presented. AS discussed in 
КООП TIE.B.2 an approximate check on the craft's roll 
stability can be nade if the SES geometry and the vertical 
lE On Of its center of gravity are known (Fig. 1). Ina 
Beet turn the principal force effecting the roll behavior 
arises from the outboard side of the starboard sidewall 
(celativ2 high pressure on the structure there due to wave 
buildup) while there is only a small force (due to snall 
pressure change) on the inboard side of the port sidewall 
Semeraputing a roll moment only if there exists a roll 
angle. Now both vertical forces are directed upward. This 


approach and the resulting forces are shown in Figure 7. 


à 1 








FZHDRP(P) FYH(S) 


Figure 6 - ACTING DEADRISH FORCES (TWO SIDEWALLS) 







FZHDRP(S) 






H(P) 
τ FYH(S) 


lu cure 7 — ACTING DEADRISE FORCES (ONE SIDEWALL) 


Bore < Since the vertical force contributed by the inboard 
Sidewall is rather small, as shown above, it can be 
E C Ceed in the roll stability check (Sect. III.B.2). 
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To investigate the effect of this interpretation of the 
acting forces two statements controlling the computation of 
deadrise forces via the sign of the rudder angle RUDANG 
(port rudder is positive) and the sidewali under 
investigation (PM1=-1 for port side) have been added to the 


computation procedure for the vertical projected force 


DDRANG=0.0 

IF (DS.GT.0.5833) DDRANG=(DS-0.5833) *0.0629 
DRANG=1.021+DDRANG-PM1* PHI 

CTN DR=COTAN (DRANG) 

RUD SIG=SIGN (1.,RUDANG) 

IF(RUDSIG.NE.PM1) CTNDR-PMT1*TAN (PHI) 

FZ OLD (J) =FZH (J) 

PZHDRP (J)zPM1*FYH(J)*CTNDR*PROMOT 

an (3) =FZH (J) +PZEDRP (J) 


The steady state values for simulated turn maneuvers using 
thrust mapping and center of gravity deadrise angle are 
shown below considering two cross-flow drag coefficients and 


both lateral added mass coefficients Uer 


TABLE Xi 


EE m m mam ae qm q m q Q s x= m <a m `V 


πμ... Ís m ae s < = шшр шшш < J m == 


Rudder angle Roll angle (deg) Speed (kn) 
(deg) simulation Ref. 6 simulation Ref. 6 
CDSW=1.28 
5 ο 5710. 52) ο ο. ΙΙΙ; 5) 9.8 
5 ντ υυ 0028 ПОСТОЕН) ο. 
12 12227012 33) Oeo τσ πας 8) Vee 2Z 
CDSW=1.16 
5 O24 10025 1) 0.09 19.3 (19:5) 1923 
9 Оа (0. 99) 08285 τ 15} 19,5 
l2 50 O 91 ino c9 9) 15:2 
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Comparing the steady state values for simulated runs in 
maples X and XI , the effect of using only the outward 
sidewall for the deadrise force Computation can be observed 
to give less favorable (little larger) steady state roll 
angles but more favorable (larger) steady state velocities. 
This is exactly what has been expected since, as shown in 
Figures 6 and 7, if deadrise force contributions from both 
sidewalls are used, the deadrise force from the port 
sidewall will reinforce the one fron the starboard sidewall 
toward a less outward roll moment (Fig. 6) while applying 
the philosophy depicted in Fig. the vertical force from 


7 
the port sidewall counteracts the starboard deadrise force. 


Regarding the roll angle responses shown in Appendix A 
for 12 degrees rudder angle with CDSW=1.28 (plots 41 and 43) 
and CDSW=1.16 (plots 45 and 47) and comparing these with the 
corresponding ones from the previous section (plot 33, 35 
and 37, 39) the change in the vertical force computation 
considering the outward sidewall only is seen to result in 
BENE 8 versus 10 cycles of transient oscillations with less 


amplitude. Using Ze (plots 43 and 47) has a auite 
Significant damping effect and makes the roll angle response 


Ecc resemble a step, which is in accordance with the 
experience gained by Ref. 7. Also, comparing these plots 
EN Plots 35 and 39, the negative roll angle occurring at 
the moment when the rudder has been introduced does not show 


up any more. 


44 





EN VEZRITTCAL LOCATION OF CENTER OF GRAVITY 


In Section III.B.2 it was discussed how to check for the 
K (C Craftis roll bahavior in a turn provided the SES 
geometry and the location of the center of gravity are 
known. This section will be concerned with the effect of 
different vertical locations of the center of gravity on the 


Emulated XR-3 roll behavior 


From experimental aeasurement the location of the center 
of gravity of the XR-3 has been determined to be 10.05 ft 
forward of the transom and Leo and Boncal established the 
vertical location at 2.54 ft above the keel on the 
memgetudinal center line of the craft. As modifications to 
the testcraft were introduced (e.g. engines and seals) the 
horizontai location of the center of gravity was again 
established experimentally. No such changes were made on the 
vertical location. In order to investigate the sensitivity 
of the steady state roll angle to vertical locations of the 
Genter of gravity, reductions in height (in accordance with 
the stability check geometry depicted in Figure 1) in 
increments of 0.1 ft were entered into the Loads and Motion 


К Сагап. 


Simulation runs with constant thrust and initial speed 
of ЕШ Кп introducing a 15 degrees port rudúer angle, 
deadrise force contributions from both sidewalls computed at 
the center of. gravity, the sidewall drag coefficient 
CDSW-1.28 and the lateral added mass coefficient being 0.4 
Were executed and the results are shown in Table XII. From 
the listed data it follows that with a lower vertical 
location of the center of gravity a tendency toward an 


inward roll angle (here actually a smaller outward roll 
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angle) can be achieved. As shown, a decrease in vertical 
distance by 0.1 ft resulted in about 0.1 degree decrease for 
both peak and steady state roll angle (almost linearly 
related). From the corresponding plots as indicated in 
Table XII a roll damping effect due to the relocation of the 
center of gravity can hardly be recognized. The number of 
cycles of transient oscillations and the steady state speed 


Mere about the same for all runs listed in Table XII. 


— = 
νο... s — < - m = s s = = == “= С _ шш === ш = = = = = = = 


w “m a Qm se πο s a =c FP" < ma sm = ош 


vertical Location Of CG experimental 
2.54 2.44 2.34 2.24 testcraft data 
steady state 
roll angie 1555 Is 52 1,22 1316 
peak 
roll angle IMA 5 
peak 
Boll rate ОЛ ΤΥ Ι8 1.9] 
y Tes of 
transients 8 8 8 8 
refer to 
plots Πο ο ο ο. σος 53-54 
steady state 
speed ο ος ος. 19.3 15.7 


Table XIII shows tne effect of changing the vertical 
location of the center of gravity for simulation runs at 
20 kn with 12 degrees port rudder angle and thrust mapping 
applied, deadrise force computation for both sidewalls at 
the center of gravity and the coefficients being 0.4 for the 
lateral added mass and CDSW=1.28 for the sidewall drag. 


Ee conditions were similar to those used in Table Z. 
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From Table XIII the effect of reducing the vertical location 
of the center of gravity in 0.1 ft increments over a 
reasonable range can be observed to also result in linear 
reductions in maximum and steady state roll angle for the 
case if thrust mapping is used. Nonlinear reductions appear 
for the peak roll rate. The number of cycles during the 
transient period (ten) is again not effected by the vertical 
relocation of τε ceentet бї gravity. Regarding the 
corresponding plots a very slight damping effect can be 
observed for reduced vertical locations. The results given 


in Table XIII are graphically represented in Figure 8 . 


TASLE XIIL 


— Düne s “s s À сс <= mmm mu pee с mm “s s m Qm mm Q = s ES О. 


vVercehicalf]ocaotion of ος experimental 
2554002 UE 3121272 testcraft data 
Steady state 
roll angle ο ο ο ο О, 0.77 ο З 
peak 
toll angle I2 5.25 1.10 
peak 
Boll rate ΙΤ ΠΝ 1.58 
cycles of 
transients 10 10 10 10 
refer to 
plots ES 50857595 59-60 
Steady state 
speed ΙΛ ;. 1 5.4 1.2 
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2.24 2.34 2.44 2.54 


MISS ROLL CONDITIONS VERSUS 
NACO CATIONTOF CENTER OF GRAVITY (TABLE ХІІІ) 
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ROLL LAMEING DUE TO VERTICAL WAVE GENERATION 


One part of subroutine SIDEWALL is concerned with the 
roll damping effect which is due to vertical wave generation 
Пиета a roll motion. The computation of the roll damping 
term was developed by Oceanics Inc. and added to the 
program by Menzel [Ref. 2]. The principles used in tne 
development cf this addition were not provided by Oceanics 
Inc., however, the apparent idea is to reduce the computed 
roll moment FK by some value which is dependent on the 
craft's roll rate and its draft, thereby on the vertical 
added mass A33. During the study of the effect of these 
additional terms a discrepancy was noted in the dimensions 


used in the calculations. 


The roll moment FK is calculated by the equation given 
in SDWL1950 and considering the expressions leading to this 


Statement it is found that 
О ЕК - 1ο. ft 


Zeh is the correct dimension for a moment vector. Now 


Considering the dimension of the roll damping tern 
dim [PROMO2*YSW*YSW*BC2*P/PI] 
which is subtracted from FK (SDWL2200) it is found 


PRCMO2 and PI are dimensionless 


СОГОН = ft 
dim [P] - rad/sec 
demepBc7].- din [AC2] 


= din [?C2 * length] 
= dim [F2 * cos(x) * length] 
ΙΙΙ; length] 
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MIRHO < BEF B = length] 
ου εν ο το, τος τς. 


So for the rcll damping term 
dim [YSW * YSW * P * BC2] 
EN таш Бес) WD". sece/ft^) . fts 
= lb. ft. sec 


which is not the correct dimension for a moment expression. 
Assuming that only the given terms may be involved in the 
roll damping calculation, it is possible to come up with the 
proper dimension if the roll rate P were squared but keeping 


ES sign, thus getting 
mo damping term - PROMOZ*YSW*YSW'*BCZ*P*ABS(P)/PI . 


Also the vertical added mass should be associated with some 
'velocity squared! - expression which is absent in the 
ExUuinal version of this part of the program but could be 
generated by squaring the roll rate P. The effect of the 
nodified expression will be an increase in the damping 
σὲ for roll rates P>1.0 which exist only in the initial 
phase of the turn maneuver (maximum roll rate in all 
ENENation runs was about 2.0) . But for most of the run 
length P is less than unity and the supposed modification 
will decrease the damping effect especially when approaching 
Steady state (P=0.0) . 


The PROMO2 - term provides a means to arbitrarily set a 
Ep damping factor which, as stated in Reference 2, has 
been found from experiments to be 316.0. Until now 
PROMO2=1.0 has been used in this thesis work. So the effect 
of changing it to the experimental value is investigated 
next. In the runs to be discussed turn maneuvers at 20 kn 
with a 12 degrees port rudder angle have been simulated with 
deadrise forces being computed at the center of gravity, 


using the experimental damping factor PROMO2=16.0 and the 
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coefficients CDSW=1.28 for the sidewall drag and E for 


the lateral added nass. 


EN hrust is held constant with deadrise force contributions 
from both sidewalls 
Comparing the obtained plots 61-62 with plots 21-22 
(PROMO2=1.0) only a slight damping effect on roll angle 
response can be observed (peak roll angle 1.450 versus 
1.479), but the number of transient oscillations decreases 
from eleven to nine. The peak roll rate decreased from 2.14 
to 2.10 for PROMO2=16.0 and the damping effect is more 


Peonounced in the roll rate time history. 


2. Thrust is mapped with deadrise force contributions fron 
both sidewalls 

Comparing the plots 63-64 with the corresponding ones 

(33-34) for PROMO2=1.0 , again only a slight damping effect 

can be noticed. The peak roll angle is again reduced from 

1.47 to 1.450 and the number of cycles of transient 


Beertlations is uneffected. 


If the damping term is changed to include the  P*ABS(P) 
expression (plots 65-66) there result eleven versus nine 
cycles of transient oscillations With little higher 
amplitudes (peak roll angle is again 1.479 as it was for 
PROMO2=1.0) . The roll rate time history also shows little 


larger amplitudes (peak roll rate is 2.14 versus 2.10) . 


3. Thrust is mapped with deadrise force contribution from 
the outward Sidewall only and lateral added nass 
coefficient ΕΓΩ 


Comparing these plots (67-68) for PROMO2=16.0 with plots 
43-44 where PROMO2=1.0 again only a slight damping effect 


in roll angle response can be noticed with the number of 


= | 





transient oscillations being uneffected. Using the P*ABS(P) 
term the slight damping eifect shown in plot 67 is cancelled 


ЕБ помп in plot 69. 


The above observations show that the experimentally 
determined roll damping factor PROMO2=16.0 has only a 
negligible damping effect. Considering the roll damping term 
again as it was added to the present XR-3 Loads and Motions 
Meogram in 1974, there was a 1/PI - factor which actually 
cancelled the PI - factor contained in the former expression 
for the vertical added mass A33S (Section III.A). Since the 
PI - factor already has been eliminated by correcting the 
A33S - expression, there is no further need to cancel it in 
the damping term by the 1/PI - term which therefore also 
could be eliminated. This change is expected to cause an 
increase in sidewall roll damping by a factor of PI and has 


been investigated by an additional run. 


4. Thrust is mapped with deadrise force contribution from 
the outward sidewall only, roil damping term being 
16 .O*#YSW*YSW*BC2*P | 


and lateral added mass coefficient er 


From plots 71-72 again only a slight damping effect versus 
meots 67-68 can be noticed. 


The investigation described in this section shows that 
tne roll damping due to vertical wave generation as included 
in subroutine SIDEWALL is not very effective even if the 
factor  PROMO2 in the damping term is increased from 16.0 to 
ENS. There also is a lack in matching dimensions 
(statement SDWL2200). 
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In all simulation runs executed during this study the 
XR-3 Loads and Motions Program as given in Ref. 2 has been 
used except for the modifications mentioned in this work. 
This simulation program contains Forbes! version [Ref. 9] 
of 'SUBROUTINE PROP! which included some revisions of the 
original version given by Leo and Boncal [Ref. 8]. Since in 
ШЕЕ Vil (Section ІІІ.В.4) the roll moment contributed by 
the propulsion system (FKP) is rather small (0.1 to 0.2) an 
additional run using Leo and Boncal's version has been 
executed under the conditions as stated for TABLE VII 
(center of gravity deadrise angle). The steady state values 
are shown below (in parantheses are the values obtained 


using Forbes' version) : 


Batch angle 0.44 (0.55) draft Sod 15291) 
roll angle i25 7101555 3) speed ПЭ 391195 37) 
ETSW 2909. 4 (821.9) ΠΗ. (+ 3227) 
PYAED - 21.3(- 20.0) R*V*AM 888.0 ( 809.2) 
FYP 97.0( 0.0) 

FKBS -3408.0(-336.4) FKSS - 3.0(- 2.7) 
FKSW ο μου 212.1) E RUND 153003 (239,9) 
FKAED - 63.5(- 56.0) FKP -305.2( 0.14) 


ABPB*PHI*Z 299.5( 273.5) 


FYD -729. 2(-651. υ) РҮН (P/S) -43/-137(-u5/-126) 
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Regarding these values drastic differences are found for 
FYP and FYRUD as well as for FKP and FKRUD, though the 
steady state values for pitch and roll angle, draft and 
speed are close to those obtained using Forbes! version. 
The pitch and roll angle responses are given in plots 25-28. 
Comparing these with the corresponding ones using Forbes! 
version (plots 13-16) it is seen that Leo and Boncal's 
version generates more oscillatory roll and pitch angle 
responses. Therefore a review of subroutine PROP especially 
the computation of the respective forces and moments has 


Been Carried out and are discussed below. 


The review resulted in force and moment equations which 
are identical to those given by Leo and Boncal (Ref. 8]. 
Ancther point that supports the supposition to use this 
version of subroutine PROP is the fact that it provides a 
reasonable roll moment due to the propulsion system (Forbes! 
Me—rslOn provides almost zero roll moment). In order to test 
me efect of changing subroutine PROP to its original forma, 
simulation runs 1. and 2. from the preceding section have 
been repeated and their time historiss for the roll behavior 
Μπ shown in plots 73-74 and 75-76, respectively. Comraring 
these with the corresponding ones using Forbes! version of 
Subroutine PROP, plots 61-62 and 63-64, it is observed that 
the steady state and peak roll angles are larger by about 
7% and the number of transien*t oscillations is unefrected 


by this change. 


Next a review of subroutine RUDDER has been carried out 
in order to investigate the cause of the change in sign for 
FYRUD and FKRUD depending on the propulsion subroutine used 
(Forbes! or Leo and Boncal's version) although identical 
RUDDER-subroutines have been used in both cases. The 


important equations to be considered are : 
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DSR=2+ZS-XR*THETA 

ENDFAC-1. * DSR/(DSR*RSPAN) 

VH=V + XR*R - ZR*P 

EFFANG=RUDANG - ENDFAC*VH/U 
FY-RHO*U*U*RAREA*ENDFAC*RCLB*EFPFANG 
° 


КЕ -AR*FY 


The above equations compute 

- draft at the rudder location 

- endfactor depending on the rudder length below the 
craft's keel 

- sideslip velocity 

- effective rudder angle depending on the craft's forward 
and sideslip velocities 

= lateral force depending on rudder area and effective 
rudder angle 


- roll moment due to rudder . 


The sign of the roll moment FK is dependent on the sign 
of the lateral force FY which depends on the sign of the 
effective rudder angle EFFANG. This in turn depends of 
course on the introduced rudder angle, but the magnitude of 
the second term of this equation  (sideslip angle) is 
responsible for the actual Sign of EFFANG. VH/U represents 
the tangent of the sideslip angle which - applying small 
mmc approximation - is subtracted from the introduced 
rudder angle (see Figure 9). Since ENDFAC=1.73 and U was 
the same tor both runs, a difference in VH which was larger 
in Leo and Boncal's version by about 50% caused the 
negative sign (versus a positive sign using Forbes! version) 
for the effective rudder angle EFFANG and thereby РҮ. So 
the only way to introduce a change toward an inward roll 
behavior here (FKRUD was positive) is to first check the 
input value for RSPAN and second to eventually modify the 


craft's engines to have larger immersion into the water. 


55 





With regard to the first item, if RSPAN would increase, 
ENDFAC would decrease as would the magnitudes of  EFFANG 


(negative), FY (negative) and FK (positive). 


U 
SSA 
ESSA 
VH 
| 
| 
| 
| 
4 
EFFANG k- RUDANG 
| 
| 
K— ESSA 
| 
SSÀ - VH/U Sideslip angle 
ESSA=ENDFAC*VH/U effective Sideslio angle 
RUDANG rudder angle 
EFFANG effective rudder angle 


Giver ΕΡΕ ANGLE WHR A TURN MANEUVER 


ti 


Figure JE I 


Considering next the fact that the engines of the XR-3 
had been replaced in 1976 a recent measurement of the rudder 
and propeller location yielded the following data (original 


input values shown in parantheses) 


КРО = -1.50 (-1.275) 
ZPO = -1.333 (-0.608) 
XRO = -1.083 (-1.125) 


PeIPAN= 1.333 ( 1.21) 
КОЕ A= 1.82 ( 0.68 ) 
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Two of the runs discussed in the preceding section have 
been repeated using the original roll damping term and the 


new propulsion and rudder input parameters. 


1. Thrust is mapped with deadrise force contributions from 
both sidewalls with Tae а 


memparing plots 77-78 with the corresponding ones for the 
former rudder and propulsion parameters (plots 63-64) the 
expected effect due to a decrease in FKRUD can be observed 
resulting in smaller peak roll angle (1.349 versus 1.459) 
and a Significantly decreased steady state roll angle (0.659 
versus 0.990 at t=50 sec). The transient period of the roll 
angle response in plot 77 shows the same number of transient 
oscillations as plot 63 but the response did not reach 
Steady state even after 45 sec . The steady state roll 


angle could be estimated to be about 0.3 to 0.4 degrees. 


2. Same conditions as in the previous run but with deadrise 
force contribution from the outward sidewall only and 
С z0.8 ° 
h 


Bearding plots 79-80 versus plots 71-72 a similar effect as 
mentioned before is observed. Now the roll angle response 
Shows a nore oscillatory transient period. The peak roll 
angle reduced from 1.46 to 1.36 degrees while the roll angle 
1-0 5ο changed from 1.34 to 0.97 degrees. For this run 
the steady state roll angle could be estimated to be about 
0.8 degrees since the roll angle response again did not 
reach steady state for the run using the recently determined 
rudder and propulsion parameters. The initial small 
negative roll angle which appeared in the previous run 
Bot 77) at the moment when the rudder motion is 
introduced, does not show up any more. fher overall roll 
response shown in plot 79 does not fall off as steeply 
(compared to plot 77) i£, in a turn maneuver, deadrise force 


contribution from the outward sidewall only is considered. 





Fron these runs it is obvious that the rudder and 
propulsion parameters of the new engines result in a more 
favorable XR-3 roll behavior with a significantly less 
steady state outward roll angle. At the time of this 
Sensitivity study Reference 7 pointed out that there 
actually is no roll angle indicated by the measuring devices 
Иа роге turn at 20 knots and 12 or 15 degrees rudder 
angle. Keeping the device sensitivity of +0.5 degrees in 
mind, it can be concluded that the steady state roll angle 
provided by the present XR-3 Loads and Motions Program 
including the before mentioned program modifications by the 
author, is quite well in agreement with the actual craft 
data. Further efforts should be undertaken to obtain a 


better damping effect during the transient period. 


Remark : 

Although the new engines extend about 0.5 ft deeper into the 
water than the old engines, the new input value for RSPAN is 
close to the former value. In Reference 1 the denominator of 
the second term in the equation for ENDFAC is defined to be 
the distance from the water surface to the bottom tip of the 
rudder. The equation for ENDFAC used in the present program 
accounts for changes in draft DSR in both the numerator and 
the denominator. The RSPAN-term shouid then be defined to 
be the distance from the craft's keel to the bottom tip of 
the rudder. The author suspects that formerly there has been 
Chosen too large a value for RSPAN in the case of the old 
engines (RSPAN=1.21 ft, probably from water surface to 
bottom tip of the rudder). Since the bottom tip of the 
rudder and the center of tne propeller were located at 
almost the Same distance from the craft's keel the 
magnitudes of ZPO and RSPAN should be about the same which 
is not true for the former chosen values. Fron these 
observations the author concludes that until now a too 
optimistic roll behavior resulted due to a too large value 
КОТ RSPAN. 
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The preceding sections of this chapter reflect the 
results of an investigation on how the XR-3 roll behavior in 
turn maneuvers at 20 kn simulated with the Loads and Motions 
Program is affected by certain changes in parameters and 
forces. The results of this sensitivity study are 


sunmarized as follows : 


x Added mass effect 


A reformulation or the vertical added mass computation which 
reduced the A33s-value by a factor of 1/PI had only a little 
effect in changing steady state draft, pitch angle and 
required thrust for constant speed in straight runs. In 
turn maneuvers this change resulted in less damped responses 
for pitch and roll angle and affected the steady state 
K for pitch angle by +8.3 % „ roll angle by -1.0 % and 
Hurt by +1.4 % . Considering the roll behavior, the effect 
Ewreducing the A33s-value by a factor of 1/PI in the 


vertical acded mass computation was negligible. 


A change in the lateral added mass coefficient о το 


0.8 had the effect of reducing the amplitudes of the 
transient oscillations and the steady state roll angle by 
0.1 degree toward the experimentally measured value. Since 
in a turn maneuver with 15 degrees rudder angle an unstable 
craft behavior showed up which had not been experienced in 


practice, the value of lee is not realistic if deadrise 
EE contributions from both sidewalls are considered. ΤΕ 


deadrise force contribution from the outward sidewall only 
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is used in the simulation program the increase in lateral 
added mass coefficient causes an increase in the steady 
state roll angle by about 0.1 degree. But since now the 
roll angle response is quite significantly damped and 
resembles a step as it does in practice, the choice of the 
larger coefficient (0.8) seems to be appropriate for the 
case of the deadrise force contributed from the outward 


sidewall only. 


* Deadrise angle effect 


The deadrise angles over the length of the curved XR-3 
Sidewalls are not uniform. Changing the deadrise angle which 
is used in the deadrise force computation from the transom 
(78.89) to the center of gravity (58.59) resulted in more 
damped responses for pitch and roll motion Dn τη 
maneuvers. For all chosen rudder angles generating port 
turns at 20 kn the steady state roll angle changed favorably 
ру 23 το 38 3 toward the corresponding roll angle measured 


experimentally and documented in Ref. 6. 


κ Cross-flow drag coefficient 


алла the cross-flow draq coefficient  (CDSW-1.28 for a 
long flat plate) to lower values, e.g. CDSW-1.00 , the 
testcraft's simulated roll behavior could be forced to match 
the neasured steady State values. Bit. Since this 
coefficient has been chosen arbitrarily and cannot ое 
declared to be appropriate for the actual XR-3 sidewalls 
(CDSW=1.0 means that there is no effective drag coefficient) 
this CDSW-value has not been considered any longer. Another 
πι «οσεθά cross-flow drag coefficient CDSW=1.16 [Ref. 3] 
resulted in about 5 4 reduction in steady state roll angle 
and a shorter transient period (about 20 %) with smaller 
amplitudes of oscillation than did the runs using CDSW=1.28. 





* Thrust mapping 

Performing thrust а eener Simulation runs using the 
data provided in Ref. 6 and assuming that both engines lose 
feat amounts of thrust in turn maneuvers, the simulated 
steady state roll angles could be shown to agree better with 
the measured angles by about another 10 % but the roll angle 
response did not quite reach steady state even after 45 


seconds. 


* Deadrise force from outward sidewall only 


The simulation runs executed with the deadrise force 
computed for the outward sidewall only and thrust mapping as 


done before, showed less agreement for both Den (0.9) 
with the corresponding measured steady state roll angles 


than did the runs considering deadrise forces from both 


sidewalls. But for URS the roll angle response was 
significantly damped and looked more like a step, as was 


experienced in practice. 


vertical location or the center of gravity 


The given value for the vertical location of the center of 
ο στ (2.54 ft above the keel) might not be quite correct 
Esto the difficulties in experimentally determining it. 
M@emesignificant effect of the vertical center of gravity 
location on the XR-3 roll behavior is obvious from Figures 1 
ES and the discussion in Section III.r. Since it controls 
Euxcecdongiy the XR-3 roll behavior it is important that 
future investigators obtain an accurate measurement of the 


Meee Location before taking testcrart data. 





* Subroutine PROP 


A review of the propulsion subroutine revealed that the 
force and moment calculations as given by Leo and Boncal 
(Ref. 8) should be used instead of the version stated by 
Forbes (Ref. 9). This change results ina 7 % larger outward 
steady state roll angle but it provides a reasonable roll 


moment due to the propulsion system. 


* Roll damping due to vertical wave generation 


τοι] damping calculation included in subroutine SIDEWALL 
has only a negligible effect on the craft's net roll 
behavior over the range of the damping factor 
1.0<PROMO2<50.0 . Since there also is a discrepancy in the 
dimension of the roll damping term the source of this part 


of the program should be investigated. 


* Suggested areas for further investigation 


During the studies undertaken with the XR-3 Loads and Motion 
Program the author found the following points either in the 
input data from the testcrait or in the program that require 


further investigation : 


a Rudder/thrust angle 

As learned from Ref. 7 the rudder/tnrust motion is directly 
introduced to the port engine while the starboard engine 
follows this via a metal bar connecting both engines. The 
EErbDoard rudder and thrust vector may be off by 2 or 3 
degrees to either side, depending on the direction of the 
Bw@eader/thrust action. Thereby the roll behavior Will 
certainly be effected. Since rudder/fthrust mapping is 
available in the Loads and Motion Program, it is important 
that accurate data be collected for each engine during a 
Bost гип. 
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а XR-3 weight 
During each test run there is a possibility that the weight 
is not exact since, as learned from Ref. 7, the testcraft 
takes on water in the seals and thereby its weight increases 
causing larger draft which also effects the roll behavior. 
This additional weight can be substantial and attempts 
should be made to measure or estimate the additional weight 


to be used in the Loads and Motion Program. 

s Cross-flow drag coefficient 
Further efforts should be undertaken to establish the proper 
value of the cross-flow drag coefficient since CDSW=1.28 
Might not be quite appropriate for the actual shape of the 
XR-3 sidewalls. 

s” Sideslip velocity 
Those parts of the Loads and Motions Program where terms 
contributing to the lateral velocity V are computed should 
be reviewed since a large value of V may result in a change 
in the effective rudder angle from positive to negative. 

In this sensitivity study it has been shown that the 
most significant effect on the simulated XR-3 zoll behavior 
are contributed by 

* deadrise angle 

* vertical location or the center of gravity 

* rudder and propeller location . 

Less significant effects are contributed by 


* thrust mapping 


EMCTOSS-rlow drag coefficient . 
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For further studies concerned with the XR-3 Loads and 
Motions Program for improved roll behavior representation it 


is recommended to use 


- the reformulated equation for the vertical added nass 
233 


- deadrise angle at the center or gravity location 


- deadrise force contribution from the outward sidewall 


only 


- subroutine PROP as given by Leo and Boncal and in 


Appendix C of this thesis 


- thrust and rudder mapping. 
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APPENDIX A 


PLOTS 


In the table given below are listed all plots with the 
respective parameters used in the corresponding simulation 


run. The table contains the following abbreviations : 


CODEF ecoefficiert CG седев отат Еу 

CDSW sidewall drag e lateral added mass coef. 
chy Chapter IV 

DRANG deadrise angle DF16 PROMO2=16.0 

PA pitch angle PR pitch rate 

RÀ roll angle RR roll rate 

REFTAB refer to table RUDANG rudder angle 

τα Transom 

= B Leo and Boncal!s version of subroutine ‘PROP! used 
EEXNP| or P*FI P replaced by ... in SDWL2200 (Section III.G) 


EN vertical location of CG is 2.585 - 0.1 (ft) 


Per REF RESP RUD DR СОЕТ CORF THRUST SBHL- REMARK 


NO TAB ANG ANG E CDSW MAP FORCES 

1 1 RA 35 TR 041027228 no both no 

2 1 РА 35 TR Оз 1.28 no both u 
B 1 RA 35 TR DT ρα. no both with 
4 1 FA 35 TR 0.4 1.28 no both E. 
5 8 RÀ 15 ER Or 1.28 no both 

6 4 FR > TUR 0.4 1.28 no both old 
7 4 PA 15 TR 0.4 1.28 no both A33s 
8 4 PR 15 TR Ong 1.28 no both 





ШОТ REF RESP RUD DR Gost ОЕР THRUSI Οκ REMARK 


NO TAB ANG ANG τ CDSW MAP FORCES 
I 

9 4,6 RA 15 TR 1228 no both 

10 4,6 RR 5 TR 0528 по both 

11 4,6 PA 15 TR 1228 no both 

12 4,6 PR 15 TR 1228 no both 

13 6 RA 15 CG 1228 no both 

14 6 RR 15 CG pug. 1.28 no both 

15 6 PA 15 CG τν ο. no Both 

16 6 PR 15 CG 1,28 no both 

17 6 RÀ 15 CG 28 no both 

18 6 RR 15 CG 28 no both 

19 6 RA 12 TR 1.23 no both 
20 6 RR 12 TR 1.28 no both 

21 6 RÀ 12 CG 1,28 no both 
p 6 RR 12 CG 1228 no both 

23 6 RÀ 12 ςς 1528 no born 
24 6 RR 12 C 1:28 no both 

25 chy RA 15 CG 025 128 no both L=B 
26 chy ER 15 CG 1.28 no Doth Lep 
27 chu PA 15 CG 1.28 no both LB 
29 chy ER 15 CG 1228 no both EB 
29 8 RA 15 CG 024 1.28 no both 

30 8 RR 15 CG 28 no both 

31 8 RA 15 CG 1. To no both 
32 8 RR 15 CG 1.16 no both 

33 10 RA 12 CG 1.253 yes both 

34 10 RR 12 CG 1228 yes botn 

35 10 RA 12 СС T28 yes both 
36 10 RR 12 ςς 0 5 1.28 yes both 

37 10 RA 12 CG pou 7216 yes both 

38 10 RR 12 CG tato ves both 

39 10 RA 12 СС 120 yes born 

40 10 RR 12 CG 15156 yes both 
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EDI REF RESP RUD DR СОСЕ СОКЕ THRUST —SDNL- REMARRK 
NO TAB ANG ANG Er CDSW MAP FORCES 
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APPENDIX B 


ОО Επ λος ο THE SIMULATION PROGRAY 


— During the course of studies undertaken with the XR-3 
Loads and Motion Program some statements have been added or 
changed in order to improve the work with the simulation 


program. 


* MAIN program 
The statements 


COMMON /AIR/ PINF,RHOINF,GAM,IQUIT (MAINO050) 
IF (MYTIME(DUM) .LT.IQUIT) GO TO 12 (MAINO66 1) 


have been modified and added to the MAIN-program. The 
purpose of the second statement is to start the writing of 
tne output values desired if computer time has reached the 
demanded time  (jobcard) reduced by IQUIT which is the 
expected amount of time (in 10 * seconds) required to print 
the output. The user may input the desired value of IQUIT as 
the fourth parameter (110-Format, columns 36-45) on the 107 
КО а сага. Default value for IQUIT is 600000 (60 seconds). 


* Subroutine INCON 
The added statements are 

IQUIT=600000 (INCNO088 1) 
to set the default value and 


READ (5,2081) IQUIT (INCN1455) 
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to read the set value for IQUIT thereby overriding the 
default value. 
* Subroutine RHS 


In Menzel's version [Ref. 2] of the XR-3 Loads and Motions 
Program the entered craft's velocity UO (kn) is changed in 
Subroutine INCON to U (ft/sec) by 


U = UO*1.6889 (INCN4580) 


mma in Subroutine RHS transformed back to units in knots 
(RHS 1850) by 


VEE = 0.5925*U . (RHS 1850) 
From these two equations follows 
MS = 0.5925*00*1.6839 


and for U0=20 kn there results VEL - 20.0135 kn. But since 
Zens desired that VELsUO the statement has been corrected 


το 
VEL = 0/1.6889 . 


Statements MAINO840, RHS 1180 have been changed accordingly. 


* Subroutine INTGRL 
ements subrcutine the statement 
MALL COLFIL (INT 1021) 


has been added in order to provide the output values already 
calculated to the user if the minimum time interval allowed 
in the integration process is undergone. Before this 


addition only the warning message 
FOU CTA TIME LESS THAN 1.0E-6 - - JOB STOPS! 


appeared and the output values have been lost. 
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APPENDIX C 


XR-3 LOADS AND MOTIONS PROGRAM 
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